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Abctract 


From  August  1989  until  August  1990,  a  simple  acoustic  telemetry  sys¬ 
tem  was  used  for  obtaining  real-time  data  from  5  Inverted  Echo  Sounders 
(lESs)  deployed  in  the  SYNOP  inlet  array  in  the  Gulf  Stream  east  of  Cape 
Hatteras.  Every  24  hours,  each  lES  calculated  a  representative  travel  time 
from  a  set  of  48  measurements  (r),  and  telemetered  that  value  to  a  lis¬ 
tening  station  on  Bermuda.  From  the  received  data,  a  daily  time  series 
of  the  depth  of  the  12°C  isotherm  (our  proxy  for  main  thermocline  depth) 
over  each  lES  was  calculated.  The  position  of  the  Gulf  Stream  North  Wall 
through  the  lES  array  was  calculated  on  a  daily  basis  from  the  thermocline 
depth  information  at  each  lES  site. 

The  telemetry  system  is  based  on  encoding  data  as  a  time  delayed  broad¬ 
cast  acoustic  signal:  the  delay  of  the  time  of  broadcast  of  the  signal,  with 
respect  to  a  reference  time,  is  proportional  to  the  data  value.  The  changes 
in  delay  time,  from  one  broadcast  signal  to  the  next,  are  recorded  at  a 
remote  receiving  station. 

The  lESs  were  recovered  in  .August  1990,  with  the  exception  of  the  one 
at  site  B2.  The  telemetered  data  from  the  lES  at  site  B2  was,  however, 
received  at  Bermuda.  The  RMS  agreement  between  thermocline  depths,  as 
calculated  from  the  data  on  tape  from  the  recovered  lESs  and  as  calculated 
from  the  received  telemetry  data,  is  20  m.  This  compares  favorably  with 
the  19  m  uncertainty  in  calibrating  the  rs  as  a  measure  of  the  thermocline 
depth.  The  RMS  agreement  between  the  position  of  the  Gulf  Stream  path 
through  the  lESs  as  calculated  from  the  tape  data  and  the  telemetry  data 
is  5  km. 

This  telemetry  system  is  not  lES  specific.  It  could  be  used  with  other 
appropriately  modified  oceanographic  instruments,  such  as  current  meters 
and  pressure  sensors. 


Statement  A  per  telecon 
Dr.  David  Evans  ONR/Code  1122 
Arlington,  VA  22217-5000 
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1  Introduction 

From  August  1989  until  August  1990,  a  system  to  acoustically  telemeter  data  from  five  modified 
Inverted  Echo  Sounders  (lESs)  moored  in  the  SYNOP  Inlet  Array  was  in  operation.  Figure  1  shows 
a  map  of  the  lES  sites.  All  but  one  of  the  lESs  were  recovered  during  R/V  Endeavor  cruise  EN216 
in  August  1990.  Although  the  lES  at  site  B2  was  lost,  its  data  was  received  at  Bermuda  throughout 
the  year-long  telemetry  experiment.  Regretfully,  the  lES  at  site  C2  did  not  perform  satisfactorily 
so  neither  the  telemetered  nor  the  recorded  data  were  usable. 

The  telemetry  scheme  is  based  on  the  simple  idea  of  encoding  the  lES  travel  time,  the  time  (r) 
taken  for  an  acoustic  signal  to  travel  to  the  ocean  surface  and  back  to  the  bottom-moored  lES,  as  a 
time  delayed  telemetered  signal.  That  is.  the  amount  of  delay  time  relative  to  some  reference  time, 
is  directly  proportional  to  r. 

In  order  to  implement  the  telemetry  scheme  three  basic  requirements  are  needed.  The  first  is  an 
lES  that  can  internally  process  the  rs  to  encode  them  as  a  delay  time  and  then  transmit  the  encoded 
signal.  The  second  is  a  mid-sound-channel  relay  station  to  receive  the  telemetered  signals  from  the 
lESs  and  retransmit  them  to  the  final  requirement,  the  distant  shore-based  receiving  station. 

The  telemetry  system  was  relatively  easy  to  construct  since  the  basic  components  already  existed 
and  required  only  relatively  minor  modifications.  The  lESs  are  already  microprocessor  controlled 
and  are  available  with  hardware  options  (described  below)  which  meet  the  systems  needs.  The 
intermediate  transponders/receivers  were  constructed  from  SOFAR  sound  sources  modified  to  receive 
as  well  as  transmit  signals.  The  base  station  was  also  already  operational  as  a  listening  station  for 
the  SOFAR  sound  sources  used  in  Tom  Rossby’s  Gulf  Stream  RAFOS  float  program  and  needed  only 
a  receiver  and  microprocessor  that  could  detect,  decode,  and  store  telemetry  signals  and  forward 
them  through  a  modem  upon  demand. 

This  report  gives  a  complete  description  of  the  telemetry  system  and  its  successful  application 
as  a  source  of  real-time  Zij  measurements  from  four  lESs  moored  edong  a  line  perpendicular  to  the 
mean  Gulf  Stream  path.  Hence,  the  telemetry  system  was  a  source  of  real-time  measurements  of  the 
GSNW  position  through  the  lES  array.  The  order  of  presentation  in  this  report  is  to  first,  describe 
the  telemetry  system;  second,  present  a  chronology  for  the  experiment;  third,  give  a  discussion  of 
the  different  methods  of  data  representation  (data  words)  for  the  internal  processor  of  the  lES,  the 
data  written  to  tape,  and  the  telemetered  data;  fourth,  describe  the  data  processing  for  both  taped 
and  telemetered  data;  fifth,  give  a  results  section;  sixth,  summarize  the  problems  encountered  in 
the  experiment  and  discuss  how  they  may  be  avoided  in  future  experiments;  and  finally,  present  a 
summary  and  conclusions  section. 


Figure  1;  Locution  and  telemetry  arrangement  of  the  five  lESs  in  the  SYNOP  inlet  array  off  of  Cape 
Hatteraa.  The  inset  is  a  schematic  of  the  data  relay  system. 


2  Telemetry  System 


The  telemetry  ssrstem  eoftifted  of  four  parts:  the  modified  lESs,  the  overhead  and  central  transpon¬ 
ders,  the  receiving  statioii  at  Bermuda,  and  the  acoustic  hook  up  (via  modems)  between  Bermuda 
and  the  University  of  RlMxle  Island  (URI).  Table  1  lists  the  deployment  and  recovery  information 
about  the  field  instruments  of  the  telemetry  system.  The  inset  of  Figure  1  shows  schematically  how 
the  system  works. 

The  basic  telemetry  scheme  transmitted  a  signal  from  the  lESs  to  Bermuda  via  some  intermediary 
receivers/transponders  moored  in  the  SOFAR  channel.  The  data  were  encoded  as  a  time-delayed 
signal  proportional  to  a  daily  averaged  r  . 


Table  1:  Field  instrument  deployment  and  recovery  information.  The  signals  from  the  original  lES 
at  site  C2  (lES  C2')  interfered  with  those  from  the  lESs  at  sites  B2  and  B.5.  The  lES  with  which  it 
was  replaced  (lES  C2)  did  not  function  properly. 


Instrument 

Site 

Deployment 

Date  Time  Cruise 

(Z)  (Z)  ID 

Recovery 
Date  Time 
(Z)  (Z) 

Cruise 

ID 

— 

Lat 

(N) 

Lon 

(W) 

Depth 

(m) 

lES  C2' 

08/12/89 

1717 

OC210 

10/17/89 

0721 

CH 

35  46.15 

73  33.00 

3450 

lES  C2 

10/17/89 

1015 

CH 

09/03/90 

0225 

EN216 

35  46.22 

73  32.75 

~3353 

lES  B2 

08/10/89 

0510 

OC210 

Not  Recovered 

35  37.01 

74  13.82 

2675 

lES  B3 

08/10/89 

0932 

OC210 

09/02/90 

1250 

EN216 

35  30.07 

74  03.40 

2960 

lES  B4 

08/11/89 

1859 

OC210 

09/02/90 

1.5.51 

EN216 

35  20  75 

74  50.60 

3325 

lES  B5 

08/11/89 

1515 

OC210 

09/02/90 

1835 

EN216 

35  12.13 

73  39.66 

3675 

CTl 

08/29/89 

1907 

OC210 

08/29/90 

0026 

EN216 

35  44,67 

73  35.47 

3480 

OTl 

08/30/89 

0416 

OC210 

08/28/90 

1558 

EN216 

35  32.80 

74  06.20 

2870 

OT2 

08/30/89 

1641 

OC210 

08/28/90 

2011 

EN216 

35  16.85 

73  44.10 

3500 

2.1  lES 

Inverted  Echo  Sounders  are  devices  which  emit  and  receive  acoustic  signals.  They  are  moored  one 
meter  above  the  ocean  floor  and  are  designed  to  transmit  acoustic  signals  towards  the  ocean  surface, 
receive  the  reflected  signed  and  record  the  round-trip  travel  time.  The  lES  concept  is  described  in 
Watts  and  Rossby  (1976). 

The  traditional  data  aquisition  scheme  of  the  lESs  follows  the  procedure  described  below.  Every 
half-hour  an  lES  transmits  a  burst  of  24  10-kHz  acoustic  signals  at  10  s  intervals.  The  receiver 
blocks  out  early  returns;  i.e.,  any  signals  which  may  arrive  before  a  time  interval  (preselected  by 
the  user)  that  is  less  than  the  shortest  time  interval  for  sound  to  travel  to  the  surface  and  back.  In 
a  twenty-four  hour  period,  48  sample  bursts  of  24  pings  (a  total  of  1152  rs  )  are  written  to  tape. 

In  addition  to  the  data  aquisition  and  internal  processing  scheme  described  above,  the  telemetry 
lESs  did  further  internal  processing  of  the  data  in  order  to  put  it  in  a  form  suitable  for  the  telemetry 
system.  Due  primarily  to  power  limitations  of  the  transponders,  only  one  signal  per  day  was  sent 
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from  each  lES  to  the  baM  station  at  Bermuda.  Thus,  the  lES  microprocessor  was  programmed  to  | 

calculate  one  representative  value  from  the  daily  set  of  measured  rs  .  For  each  half  hourly  burst 

of  24  signals,  the  microprocessor  stored  in  an  accumulator  the  value  ri^  which  corresponded  to  the  9 

first-quartile  point  in  the  distribution  of  the  rs  .  The  value  was  also  written  to  tape  along  with 

the  24  rs  from  that  burst.  The  daily  representative  r  was  calculated  after  each  24  hrs  as  the  mean 

of  the  48  values  of  ri,  stored  in  the  accumulator  (rf^)-  This  daily  representative  travel  time  was 

encoded  as  a  time-delayed  acoustic  signal  and  sent  to  a  relay  transponder.  9 

A  flow  chart  of  the  data  aquisition  scheme  is  shown  in  Figure  2.  The  optional  hardware  chosen 
for  the  Pacer  Systems  Inc.  Model  1665  lESs  used  in  the  deployment  included: 

•  The  high  output  driver  and  transducer  option  to  ensure  that  the  telemetered  signal  would  reach 

an  intermediate  transponder  (the  horizontal  component  of  the  signal  is  larger  with  this  option)  9 

•  The  50,000  /ifd  capacitor  option  and  high  current  capacity  lithium  batteries  were  used  to  ensure 
sufficient  power  while  transmitting  the  telemetered  signal  (50  ms  transmission  duration). 

2.2  Transponders  ^ 

The  lES  is  not  designed  to  send  acoustic  signals  over  long  horizontal  distances.  Its  high  frequency 

acoustic  signals  are  inefficient  for  long  range  propagation.  Moreover,  moored  near  the  sea  floor. 

the  lES  is  not  positioned  for  its  signals  to  take  advantage  of  the  SOFAR  channel.  As  illustrated 

in  Figure  3,  when  moored  at  a  nominal  depth  of  4500  m  the  lESs  10-kHz  acoustic  signal  has  a  ^ 

horizontal  range  of  about  20  km  before  intersecting  the  SOFAR  channel  (Chaplin,  1990).  Therefore, 

intermediate  transponders  moored  at  the  depth  of  the  deep  sound  channel,  ate  required  to  relay  the 

telemetered  signals  from  the  lESs  to  the  base  station.  Figure  4  shows  a  diagram  of  the  moorings  for 

the  intermediate  transponders.  A  transponder,  sitting  in  the  SOFAR  channel,  capable  of  receiving  ^ 

10-kHz  signals  and  then  retransmitting,  after  frequency  down-shifting  for  more  efficient  horizontal 

transmission,  can  service  several  lESs  within  a  radial  distance  of  20  km. 

The  telemetry  lES  sites  used  in  this  deployment  were  situated  such  that  one  intermediate 
transponder  could  service  each  pair  of  lESs.  The  inset  of  Figure  1  shows  a  schematic  of  the  the  ^ 

actual  deployment.  Two  types  of  intermediate  transponders  were  used.  There  was  one  modified 
SOFAR  float  (hereafter  denoted  CT  or  central  transponder)  that  directly  serviced  the  lES  at  site  C2 
and  relayed  the  data  from  all  other  lESs  to  Bermuda  via  260  Hz  signals.  The  other  two  intermediate 
transponders  (hereafter  denoted  as  OT  or  overhead  transponders)  each  relayed  data  from  two  lESs  ^ 

to  the  CT  at  780  Hz.  The  OTs  were  used  because  it  was  originally  believed  that  they  would  be 
less  expensive  than  the  CTs.  We  have  subsequently  found  that  the  cost  difference  is  minimal.  W'e 
would  probably  not  use  OTs  in  a  future  deployment  but  rather  all  CTs.  The  range  of  the  780  Hz 
OT  to  CT  link  is  reliable  to  about  200  km.  The  actual  distances  between  the  transponders  in  this  ^ 

deployment  were  less  than  60  km  (Chaplin,  1990).  At  each  stage  of  the  relay  system  the  frequency 
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Figure  2:  lES  data  aquisition  and  internal  proceasing  flow  chart.  In  a  one  day  cycle  1152  ra  are 
written  to  tape,  48  firat-quartile  points  ri,  are  determined,  and  1  is  telemetered. 


of  the  signal  is  down-shifted  while  preserving  its  time  delay  or,  equivalently,  the  data  information 

2.3  Base  Station 

The  base  station  was  located  at  the  U  S.  Navy  Underwater  Systems  Center  (NUSC)  field  station 
on  Bermuda  known  as  the  Tudor  Hill  laboratory.  As  previously  mentioned,  the  base  station  was 
originally  established  to  monitor  the  performance  of  the  moored  SOFAR  sound  sources.  For  that 
purpose,  occasional  down  time  of  the  station  was  not  a  critical  problem.  The  station  would,  at 
predetermined  times,  listen  for  the  sound  source  signals  and  record  the  time  of  arrival.  This  allowed 
any  clock  drifts  of  the  sources,  and  hence  RAFOS  float  position  errors,  to  be  corrected.  Since  the 
clock  drifts  were  essentially  linear  over  the  periods  normally  experienced  for  listening  station  down 
time,  the  listening  gaps  did  not  affect  the  quality  of  the  RAFOS  float  data. 

In  principle,  the  only  necessary  modification  to  the  base  station  for  its  use  in  the  lES  telemetry 
system  would  have  been  to  increase  its  number  of  daily  “listening  windows"  for  acoustic  signals. 
However,  we  made  several  other  changes  to  upgrade  the  system.  First,  we  added  the  capability  to 
store  and  transmit  data  over  modems  and  telephone  lines  back  to  our  offices  in  Rhode  Island.  The 
clock  that  the  old  receiver  used  had  a  drift  problem  which  had  to  be  closely  monitored.  For  this 
reason  a  completely  new  receiving  system  was  installed  at  the  station.  To  ensure  redundancy,  the 
existing  receiving  system  (RO)  was  left  intact  to  operate  in  parallel  with  the  new  one.  As  with  the 
old  receiver  clock,  the  clocks  of  the  new  receivers  were  periodically  checked  against  a  very  stable 
Cesium  clock  at  the  Naval  lab. 
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Figure  4'.  Schematic  of  Intermediate  Transponder  mooring. 


The  new  system  had  two  different  receiving  units,  a  hetrodyne  receiver  (Rl),  like  the  original 
receiving  system,  and  a  homodyne  receiver  (R2).  The  homodyne  receiver  was  installed  to  test  its 
performance  relative  to  the  hetrodyne  receiver  for  possible  use  in  future  systems.  For  each  listening 
period,  each  receiver  would  select  the  two  signals  most  correlated  with  a  reference  signal.  Unlike 
the  old  system  which  just  recorded  the  arrival  times  of  detected  signals,  the  new  receiving  system 
automatically  converted  each  arrival  time  into  the  time  delay,  in  deciseconds,  between  the  arrival 
time  and  a  reference  time.  The  lime  delays  as  well  as  the  correlation  values  of  the  two  selected 
signals  were  available  in  the  data  record.  The  new  system  had  the  capability  to  store  and  transmit 
data  via  modems.  This  made  it  possible  to  transfer  the  data  to  computers  at  URI  at  convenient 
intervals.  The  data  from  each  receiver  was  also  sent  to  a  printer  so  that  the  Navy  could  monitor 
what  was  being  sent  out  from  the  lab.  Figure  5  shows  a  block-diagram  of  the  receiving  system  on 
Bermuda. 

2.4  Signal  Timing 

Figure  6  shows  the  timing  of  the  telemetered  signal  as  it  is  relayed  from  an  lES  to  Bermuda.  AT  is 
the  0-8.5  min-delay-encoded  signal,  of  50  ms  duration,  transmitted  by  the  lES.  After  detecting  the 
signal,  the  OT  retransmits  the  signal  to  the  CT  AToc.  the  delay  time  between  transmission  from 
the  OT  until  reception  at  the  CT,  and  ATcs-  the  delay  time  between  transmission  from  the  CT 
and  reception  at  Bermuda,  are  nearly  constant  because  they  are  primarily  related  to  the  distance 
separating  the  sites.  The  time  of  arrival  of  the  signal  at  Bermuda  is  linearly  related  to  AT.  Tom 
Rossby  (personal  communication)  had  found  earlier  that  the  scatter  in  arrival  times  of  a  SOFAR 
sound  source  moored  south  of  Cape  Hatteras  was  approximately  2  s.  This  is  the  nominal  scatter  we 
expected  for  “good”  signals  received  at  Bermuda. 

The  multiple  signals  from  both  the  lES  telemetry  system  and  the  RAFOS  navigation  sound 
sources  followed  a  fixed  transmission  schedule  so  that  there  was  no  interference  between  them  at 
the  listening  station  on  Bermuda.  The  schedule  for  lES  telemetry  transmission,  OT  and  CT  listen¬ 
ing/transmission,  RAFOS  navigation  transmission,  and  base  station  listening  is  shown  in  Figure  7 

3  Chronology 

A  brief  chronology  of  the  telemetry  deployment  is  shown  in  Table  2.  Note  that  the  telemetry  signals 
from  the  original  lES  at  site  C2  were  interfering  with  those  from  sites  B2  and  B5.  Thus,  telemetered 
data  from  the  lESs  at  sites  B2  and  B5  were  not  received  until  after  the  lES  at  site  C2  was  changed. 
The  resulting  ~50  day  gaps  at  the  beginning  of  the  telemetry  records  for  lESs  B2  and  B5  are  seen 
in  Figure  15. 
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Figure  5:  Block  Diagram  of  Base  Station  on  Bermuda.  The  homodyne  receiver  was  evaluated  for 
possible  use  with  the  RAFOS  floats.  It  proved  to  be  as  effective  as  the  hetrodyne  receiver.  Note 
that  the  modem  could  only  be  connected  to  one  receiver  at  a  time. 
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Figure  6;  The  timing  of  the  delay  encoded  telemetered  signal  as  it  is  relayed  from  an  lES  to  Bermuda. 
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Figure  7:  Schedule  for  telemetry  system  and  RAFOS  navigation  reception/transmission 
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Table  2:  Chronology  of  the  Telemetry  Deployment 


Date _ 

Aug  10,  1989 


Mid.  Aug 
Aug  ? 

Oct  3,1989 
Oct  16, 1989 
Oct  17,  1989 

Nov  30.1989 
Jan  15-16.1990 

.■\ug-Sept  1990 


Comment 


The  receiver  (Rl)  was  installed  on  Bermuda  by  Skip  Carter.  The  old 
receiver  (R0),which  was  connected  to  the  ticker-tape  printer 
was  accidently  unplugged. 

lESs,  modified  to  telemeter  their  data  to  Bermuda,  were  deployed. 
Old  Ce  clock  failed  sometime  in  August. 

Receiver  RO  reattached  and  reset 
New  Ce  clock  installed 

Replaced  FES  C2  (Gerry  Chaplin  on  Cape  Henlopen).  Started 
receiving  from  lESs  B2  and  Bo. 

Jim  Eontaine  resets  the  Bermuda  Ce  clock  and  all  of  the  other  clocks. 
Reset  clock  in  receiver  Rl  (Chaplin  tests  the  homodyne  and 
hetrodyne  (receivers). 

Telemetry  lESs  pulled  out  of  water. 


4  Data  Representation:  Telemetered  And  Taped 

In  this  section  an  overview  is  given  of  the  steps  in  the  evolution  of  an  lES  internal  data  word 
into  a  data  word  written  to  tape  and  a  telemetered  data  word.  In  the  former  case  the  internal 
data  word  represents  an  "instantaneous"  r  .  while  in  the  latter  case  it  represents  a  daily  average 
T  (See  Section  2.1).  In  order  to  lay  the  ground  work  for  subsequent  sections,  which  will  include 
data  processing  descriptions  and  comparisons  between  taped  and  telemetered  data,  the  conversion 
of  both  types  of  datawords  into  thermocline  depth  will  also  be  described 

The  lES  internal  clock  is  nominally  set  to  20480  Hz.  Each  count  (lESCT)  of  its  data  word  is 
then  equivalent  to  20480“'  s  which  is  the  time  resolution  of  the  instrument. 

For  an  lES  deployed  at  5000  m  depth,  the  travel  time  measurement  is  nominally  6  9  s  or  141241 
lESCTs.  In  general,  18  bits  are  required  to  store  the  full  data  word.  However,  the  data  word 
written  to  tape  is  the  lES  data  word  truncated  to  the  13  least  significant  bits  (LSB).  This  effectively 
subtracts  off  a  common  amount  of  time  from  the  measured  rs  ,  which  is  due  to  the  depth  of  the 
instrument  not  to  changes  in  the  average  temperature  of  the  water  column.  The  conversion  between 
lESCTs  and  the  data  written  to  tape  (SCTs)  is 

lESCT  =  iVCr  +  .\CTo  ( 1 ) 

and  the  conversion  between  r  and  .\'CT  is 

^  NCT  +  SCTo 

’’  20480//r  '  ’ 

where,  NCTq  represents  a  standard  amount  of  time  for  each  lES  which  is  subtracted  off  by  truncating 
the  data  word.  The  time 

,  _  SCT 

20480//-- 


(3) 
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can  b«  thought  of  aa  a  travel  time  anomaly. 


The  travel  time  r ,  meaaured  by  the  lES,  ia  linearly  related  to  the  depth  of  the  main  thermocline, 
which  we  define  aa  the  depth  of  the  12‘C  iaotherm  (Zu  ).  The  relationahip 


,  NCT+NCTo  .  „ 

=  ”  2MS0H, 


m 


can  be  written 


where 


Zi2  =  m 


NCT 

20480WZ 


+  B' 


fl'  =  fl  +  m 


NCTo 
20480// r 


(5) 

(6) 


Thus,  when  calibrating  an  lES  using  r'  a  different  intercept  is  calculated  than  if  the  original  data 
word  T  was  available.  However,  both  calibrations  give  identical  Z12,  so  it  is  not  necessary  to  keep 
track  of  the  constant  NCTq. 

The  telemetered  data  word  manipulation  and  subsequent  encoding  is  somewhat  more  involved, 
but  it  is  still  straightforward.  Before  encoding  the  travel  time  into  a  delay  time,  the  lES  truncates 
the  data  word  from  the  full  18  bits  to  the  middle  8  bits,  keeping  only  bits  4-11  (See  Figure  8).  Thus, 
not  only  is  a  common  time  subtracted  off  of  each  r  by  throwing  out  the  7  most  significant  bits 
(MSB),  but  each  count  (TELCT)  of  the  new  data  word  is  equal  to  eight  of  the  original  counts  since 
the  3  LSBs  are  thrown  out.  The  data  word  manipulation  is  schematically  shown  in  Figure  8.  The 
conversion  between  TELCTs  and  lESCTa  is  then 


lESCT  =  8  •  TELCT  -f  lESCTo, 


(7) 


where  IESCTq  is  a  constant  number  of  counts  subtracted  off  by  truncating  the  7  MSBs. 

The  lES  constructs  a  time  delay  (At)  which  is  equal  to  TELCT*1  s,  i  e.,  with  a  granularity  of  2 
s.  Thus,  2  s  of  delay  is  allowed  for  each  8  counts  of  the  original  data  word.  It  is  important  to  note, 
however,  that  this  conversion  is  only  accurate  to  7  lESCTa  because  of  the  truncation  of  the  least 
significant  bits  of  the  original  data  word.  The  resolution  of  r  has  been  reduced  from  (20480  Hz)~'  to 
8  ■  (20480//z)~‘  because  of  the  truncation  of  the  3  LSBs  from  the  original  data  word.  Equivalently, 
the  resolution  has  been  reduced  from  1.0  m  to  7.8  m  of  Z13  (the  slope  m  in  the  above  equation  has 
been  found  to  be  -19840  m  s“‘) 

At  Bermuda,  the  signal  is  measured  in  deci-seconds  of  delay  time  of  the  received  signal.  That 
is,  each  count  of  the  data  word  (BDACT)  is  equal  to  1  deci-second  of  delay  time.  The  conversion 
between  BDACT  (received  delay  time)  at  Bermuda  and  TELCTa  is 

- -  Ida  - (BDACT -^BDACTo) 

=  10dss->-2s 

where  BDACTo  is  a  fixed  positive  or  negative  quantity  representing  delay  time  offset  which  arises 
because  the  reference  time  on  Bermuda,  for  the  telemetered  signal,  is  not  exactly  equal  to  the  sum 
of  the  reference  time  at  the  lES  and  the  signal  transit  time  from  the  lES  to  Bermuda. 


IE5CT" 


"NCT 


-TELCT 


Figure  8:  All  of  the  lES  internal  processing  is  done  with  18-bit  words.  The  data  word  is  truncated 
to  the  13  least  significant  bits  before  being  written  to  tape.  Before  being  encoded  as  a  time  delayed 
telemetered  signal  the  data  word  is  truncated  to  bits  4-11. 
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Table  3;  The  converaioiu  between  r,  the  counts  written  to  the  lES  tape  {NCT)  and  the  counU 
measured  at  Bermuda  (BDACT).  Ab  is  a  constant  for  each  lES  and  is  introduced  by  any  travel 
time  between  the  lES  and  Bermuda  of  the  telemetered  signal  which  is  not  accounted  for  in  the  time 
delay  introduced  between  when  the  signal  is  sent  from  the  lES  and  when  the  Bermuda  receiving 
station  begins  to  “listen” . 


r 

(s)10-3 

NCT 

Delay  Time 

(s) 

BDACT 

(ds) 

2 

4 

0.5 

5  -1-  No 

25 

512 

64 

640  +  No 

50 

1024 

128 

1280  -1-  No 

In  a  manner  similar  to  equation  o  for  the  data  written  to  tape,  Z\2  can  be  related  to  BDACTs 
through  the  relation 

Zn  =  m  C  BDACT  A-  (9) 

where  C  is  the  conversion  factor  between  BDACTs  and  travel  time  of  the  lES  measuring  signal, 
and  B^  is  the  modified  intercept  which  has  absorbed  both  the  truncation  of  the  MSBs  and  constant 
delay  time  offsets  {BDACTq).  The  conversion  factor  C,  is  given  by 


C  = 


_ Ids  8 _ 

lOdss-*  2s  ■20480//I’ 


(101 


or 


C  = 


Is 


51200 

Note  that  C  for  the  tape  data  is  (20480  Hz)"*. 

Table  3  lists  the  lES  data  written  to  tape  and  telemetered  data  for  a  few  sample  r  values. 
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5  Processing 


This  section  describes  the  processing  for  the  telemetered  and  taped  data.  Flow  charts  of  the  pro¬ 
cessing  steps  for  taped  and  telemetered  data  is  shown  are  Figure  9.  These  charts  assume  that  the 
lESs  have  been  calibrated,  and  in  the  case  of  the  telemetered  data,  that  they  have  been  de-drifted. 

The  details  of  the  processing  steps  are  described  in  the  following  subsections.  Note  that  the  taped 
data  consists  of  the  “traditional"  data  (a  burst  of  24  rs  every  1/2  h)  and  data  used  for  telemetry 
(the  representative  r  (riy)  calculated  by  the  lES  for  each  burst). 

5.1  Taped  Data 

The  processing  of  the  taped  data  is  fully  described  in  Fields,  et.  al.,  (1991),  so  the  description 
here  will  be  brief.  As  described  in  Section  2.1,  every  half  hour  an  lES  emits  a  burst  of  24  acoustic 
measuring  signals.  The  lES  writes  all  24  rs  to  tape  along  with  a  time  which  corresponds  to  the  # 

first  quartile  point  in  the  distribution  of  the  24  rs  .  The  first  24  data  words  (rs  )  are  used  in  the 


) 
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Taped  Data 
{ 1/2-hourly  record } 


1/2  hourly  rsoord  of  burst  of 
24  x’S  ( {NCTs) .  1/2-hr ) 


Convert  each  i/2  hourty  bursty 
of  24  X  'S  into  1  X  / 


Fill  in  missing  samples  by 
linear  interpolation  . 


1/2  hourly  record  of 
tin's  ( NCT-s  .  1/2-hr) 


'  Create  24-hr  sample  rate  by 
averaging  every  48 1,  „ 's  > 


'Convert  each  24-hr  average  x > 
into  a  delay  fime  / 


/  Convert  NCT's  into  seconds  / 

z'  N 

/  (S,  1/2-hr)  / 

^ - "T - 

1  Delay-time  record  i 
1  (  BDACTs  .  24-hr )  J 

Remove  the  tidal  effects  in 
the  1/2  hourly  record 


Oespike  the  record 


'Correct  for  seasonal  variation 
of  upper  water  column  / 


Apply  40-hr  lowpass  filter 


Convert  X 's  into  2,2 'S 


'Subsampie  data.  Reduce  frony 
1/2  hourly  to  6  hr  record.  / 


Compare  with 
telemetered 
V  record  ? 


Thermocline  depth  record 

(2,2.6-hr) 


"Comparison  1' 


Create  24-hr  sample  rate  by 
averaging  every  4  Zi2'5  / 

ihermocliriedepthrecoid\ 
(2,2. 24-hr)  H 


Telemetered  Data 
{  24-hourly  record } 


Delay- time  record 
(  BDACTs  .  24-hr ) 


/  Despike  the  record  / 


/'Correct  lor  seasonal  variation  / 
of  upper  water  column  temp.  / 


/  Convert  T 's  into  2i2's  / 


/  Fill  in  the  data  gaps  / 


/  /\ppiy  4-day  lowpass  filter  / 


Thermocline  depth  record 
(2,2. 24-hr) 


"Comparison  2 


Figure  9:  The  flow  charts  showing  the  basic  processing  of  the  taped  and  telemetered  data  are  shown 
on  the  left  and  right  panels  respectively.  Also  indicated  are  the  types  of  comparisons  made  between 
the  taped  and  telemetered  data. 
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traditional  procewing  of  the  lES  taped  data  and  the  25th  data  word  (n,  )  is  used  exclusively  for 
reconstructing  the  data  telemetered  from  the  lES. 

5.1.1  Calibration  of  r  to  Zu 

The  taped  records  are  converted  into  Ztj  records  using  Equation  5.  The  factor  m  in  the  equation  is 
known  in  the  Gulf  Stream  region  (see  Howden,  et.  al.,(1992)),  and  B'  (or  Bintiap,)  is  determined 
individually  for  each  lES.  The  calibration  curve  for  each  lES  is  determined  once  its  Binttapi  is  found 
The  following  are  the  steps  for  finding  Binttape  for  each  lES.  During  deployment  and  recovery  cruises, 
a  total  of  Ni  XBT  casts  are  taken  at  the  ith  lES  site  The  temperature  profile  of  each  XBT  cast 
is  integrated  from  250  m  to  750  m  yielding  a  quantity  Q  {Q  =  f^^Q^Tdz).  The  quantity  Z'j  is 
calculated  from  an  emperical  functional  fit  of  Zn  to  Q  from  over  6,000  XBT  casts  (from  SODC 
archives)  taken  in  the  Gulf  Stream  region.  The  .V,  Zjj  values  are  substituted  into  Equation  5  along 
with  the  .V,  coincident  r  measurements  from  the  recovered  data  tape,  yielding  S,  values  of  Bint,ape 
The  Binttape  used  in  the  calibration  curve  is  the  average  of  the  .V,  BinttapeS- 

The  motivation  of  using  Z'2,  rather  than  Z12,  in  the  calibration  process  is  that  it  eliminates 
the  high  frequency  noise  that  can  be  present  when  a  point  measurement  (Z13  )  is  made  of  the 
thermocline  depth.  Hereafter  Z*2  will  be  simply  denoted  as  Zu,  it  being  understood  that  in  the 
calibration  process  Z12  refers  to  the  calculated  value  ZJ2 

5.1.2  Traditional  Processing  of  Taped  lES  Data 

After  an  lES  has  been  recovered  and  its  data  tape  read  into  a  computer,  a  number  of  preliminary 
steps  are  taken  to  obtain  a  record  of  r  versus  time  (see  Fields,  et.,  al.,  (1991)).  For  this  discussion, 
the  important  step  is  the  fit  of  each  half-hourly  set  of  24  rs  to  a  Rayleigh  distribution.  The  mode 
of  each  distribution  is  taken  to  be  r  for  that  particular  sample.  Unlike  the  telemetered  data  record 
which  is  composed  of  daily  averaged  rs  ,  the  taped  record  has  a  half-hourly  sample  rate  and  has 
to  be  corrected  for  changes  in  travel  time  due  to  the  tides.  Next,  the  taped  record  is  passed 
through  a  program  which  removes  outliers  and  replaces  them  with  linearly  interpolated  values.  The 
“cleaned  up”  data  is  then  corrected  for  changes  in  the  travel  time  due  to  the  seasonal  variation  of  the 
temperature  of  the  upper  water  column.  A  low-pass  2nd  order  Butterworth  filter,  with  a  2-d  cutoff 
]>eriod,  is  then  applied  forwards  and  backwards  to  the  data  to  remove  high  frequency  fluctuations. 
The  record  is  then  subsampled,  reducing  it  from  a  1/2  hour  sample  rate  to  a  6  hr  sample  rate. 
Finally,  the  r  record  is  converted  to  a  Z13  record  using  Equation  5  and  the  Bints  in  the  third  column 
of  Table  4. 

For  comparison  with  the  (Z^,  24-hr)  telemetered  record  {Comparison  2),  the  fully  processed 
(Z12,  6-hr)  record  is  converted  to  a  24-hr  sample  rate  by  averaging  every  4  Zijs  .  This,  of  course,  is 
not  part  of  the  traditional  processing  but  is  simply  done  for  the  compairison. 
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5.1.3  Reconitruction  of  the  Telemetered  Data 

As  described  in  Section  2.1,  the  25th  data  word  of  each  taped  record  is  the  representative  travel  time 
(ri,  )  of  the  burst  of  24  rt  in  the  record  chosen  by  the  lES.  Internally,  the  lES  averages  every  48 
of  the  these  data  words  and  calculates  a  delay  time  proportional  to  the  averaged  value.  To  convert 
the  1/2-hourly  record  of  ri,s  into  the  information  that  was  telemetered  from  the  lES,  the  following 
steps  are  done.  First,  every  48  ri,s  are  averaged  to  produce  a  24-hr  sample  rate  (rf7,  24-hr).  The 
are  then  converted  from  NCTs  into  TELCTs  by  dividing  by  8.  since  each  TELCT  equals  8 
SCTs.  Next  a  delay  time  is  calculated  using 

Ar=rI7.20ds  (12) 

Tlie  resulting  record  (AT.  2i-hr)  can  tlien  be  directly  compared  with  the  delay  times  as  measured 
at  Bermuda  (Comparison  ]). 

5.2  Telemetered  Data 

The  basic  processing  of  the  telemetered  data  is  done  in  the  6  steps  that  are  described  in  the  following 
subsections.  Figures  10  to  13  show  the  raw  data  and  the  results  of  each  processing  step  that  give  a 
visually  distinct  change  between  the  input  and  output. 

5.2.1  lES  Clock'Drift  Correction 

Because  of  inaccuracies  in  the  setting  of  the  lES  clock  frequencies,  the  telemetry  records  had  to  be 
de-drifted.  The  clocks  used  in  the  lESs  are  accurate  to  0.01  s  d~'.  A  drift  rate  of  this  magnitude 
was  deemed  small  enough  to  yield  ihcrmoclinc  depth  estimates  within  our  stated  accuracy  of  20  m. 
For  example,  a  drift  rate  of  0.01  ds  d“*  would  yield  a  maximum  clock  drift  of  3.65  s  over  the  one 
year  deployment  of  the  instruments  which  in  turn  would  translate  into  an  apparent  drift  of  Zi2  of  14 
m.  Although  the  lES  r  measurement  has  a  resolution  of  1  m,  the  accuracy  of  an  lES  calibration  into 
Zi2  is  only  20  m  because  of  the  uncertainty  of  the  intercept  (Bint)  of  the  linear  relation  between  Z12 
and  r  .  Since  the  calibration  uncertainty  is  greater  than  the  error  introduced  by  an  actual  frequency 
drift  of  the  lES  clocks,  we  elected  to  simplify  the  telemetry  design  by  excluding  a  method  to  check  for 
clock  drifts.  Unfortunately,  we  discovered  after  the  experiment  was  underway  that  the  instrument 
used  to  set  the  lES  clock  frequencies  was  not  accurately  calibrated.  Thus  the  clock  timing  effectively 
drifted  at  a  much  higher  rate  than  anticipated.  While  these  clock  drifts  do  not  significantly  affect 
the  actual  travel  times  measured  by  the  lESs,  they  do  contaminate  the  telemetered  data.  The  basic 
information  of  the  telemetry  system  is  in  the  arrival  time  of  the  signal  relative  to  a  fixed  reference 
time.  Large  clock  drifts  would  then  shift  the  arrival  times,  which  would  in  turn  be  falsely  interpreted 
as  a  change  in  Z12  . 
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BDACT  {da) 


Z,2  Depth  (m)  Z.j  Depth  (m)  Z,2  Depth  (m) 


Figure  12:  The  top  panel  shows  the  BDACT  as  measured  at  Bermuda  for  the  lES  at  site  B4 
Successive  panels  show  the  output  of  the  processing  steps  which  are  visually  distinct. 
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While  the  lESs  were  still  moored  at  sea,  the  frequency  counter  was  recalibrated  and  crude 
estimates  of  the  clock  drifts  were  made.  This  allowed  us  to  make  preliminary  corrections  to  the 
telemetered  data  while  the  system  was  still  operational.  However,  more  accurate  drift  corrections 
could  not  be  performed  until  the  lESs  were  recovered.  The  drift  rate  was  calculated  individually 
for  each  lES  as  the  slope  of  the  best  straight  line  fit  to  the  difference  between  the  tape  records  and 
telemetered  data  records  after  outliers  had  been  windowed  out.  An  example  of  this  process  is  shown 
in  Figure  14  for  lES  B3.  Since  lES  B2  was  not  recovered  we  were  unable  to  use  this  technique  at 
that  site.  Instead,  we  simulated  the  B2  tape  data  by  averaging  the  B1  and  B3  tape  records.  Then 
we  compared  the  simulated  data  with  the  actual  telemetered  B2  data  to  obtain  the  drift  rate  of  the 
clock. 


5.2.2  Calibration  of  r  to  Zi; 


In  principle,  the  Bints  for  the  telemetered  data  flints, ,i  are  calculated  identically  to  those  for  the 
tape  records:  the  Zjo  obtained  from  XBT  casts  at  the  lES  sites  are  compared  with  the  coincident 
delay  times  (encoded  r  )  received  at  Bermuda.  However,  the  telemetry  system  was  not  operational 
when  the  lESs  were  first  deployed  and  the  XBT  casts  taken.  Thus  a  crude  calibration  was  initially 
done  using  the  XBT  casts  taken  at  deployment  and  the  first  telemetered  rs  received  at  Bermuda 
several  days  after  deployment.  Subsequently.  af.i.r  me  lESs  were  recovered,  we  were  able  to  perform 
better  calibrations  by  using  the  ♦aped  uata  which  were  calibrated  with  coincident  XBT  casts  taken 
during  both  the  deployment  and  recovery  cruises.  The  Bints  for  the  de-drifted  telemetry  records 
were  calculated  in  this  manner  as 


Bint, el 


—  19840ms”* 


(BDACT,af,e  -  B PACT, el) 

51200dss-‘ 


+  Bint, ape. 


(13) 


where  BDACT,ei  are  the  data  words  of  the  received  telemetered  data  and  BDACT,ape  are  the 
converted  taped  data  words. 

The  drift  rates  are  shown  in  Table  4,  along  with  the  Bints  for  both  the  telemetry  data  records 
and  the  tape  data.  The  table  shows  two  entries  for  the  telemetry  Bints.  The  first  entry,  the  direct 
method,  is  the  Bint  calculated  from  Equation  9  using  the  Zijs  of  the  XBT  casts  taken  during  recovery 
and  deployment  and  the  first  and  last  rs  from  the  telemetry  records.  Because  the  XBT  casts  and  rs 
were  separated  by  as  much  as  several  days,  this  calibration  is  not  as  accurate  as  that  for  the  taped 
data  for  which  coincident  XBT  casts  and  rs  are  available.  At  deployment  this  method  must  be  used 
to  obtain  an  initial  calibration  using  the  XBT  casts  taken  at  that  time.  In  a  future  deployment,  it 
would  be  desirable  to  have  the  telemetry  system  operational  so  that  the  initial  calibration  could  be 
done  with  coincident  (i  e.,  within  one  day)  XBT  casts  and  telemetered  rs  As  discussed  in  Section  3, 
the  telemetered  signals  for  lESs  B2  and  B5  were  blocked  by  the  signals  from  the  original  lES  at  site 
C2  for  the  first  ~46  days  of  deployment.  The  initial  calibrations  for  lESs  B2  and  B5  were  inferred 


BDACTS,ope  -  BDACTS,ei  (ds)  BDACTSjpp^  -  BDACTS 
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from  the  initial  calibrations  of  lESs  B3  and  B4  using  the  following  method.  At  deployment,  Zij  was 
known  at  sites  B2  and  B5  from  XBT  casts  taken  at  that  time.  Although  coincident  rs  were  not 
telemetered  from  these  two  lESs,  relatively  coincident  (i.e.,  within  several  days)  rs  were  telemetered 
from  adjacent  lESs  (B3  and  B4).  Once  the  lESs  at  sites  B2  and  B5  began  telemetering  their  data, 
the  average  and  the  rs  from  the  adjacent  lESs  were  computed.  These  differences  were  added  to  the 
coincident  rs  from  the  adjacent  sites  to  obtain  estimated  coincident  rs  from  the  lESs  at  sites  B2 
and  B5. 

The  second  column  in  Table  4  contains  the  Bints  calculated  by  comparing  the  telemetry  record 
and  the  calibrated  tape  record  as  shown  in  Equation  13.  There  is  no  entry  for  the  Bint  calculated 
by  the  comparison  method  for  lES  B2  since  the  instrument  was  not  recovered. 


Table  4:  BINTs  and  drift  rates.  Columns  2  and  3  show  the  BINTs  for  the  telemetered  data  calculated 
by  two  different  methods.  The  first  method  uses  the  coincident  rs  ,  from  the  lESs,  and  Z12  from 
XBT  casts.  The  second  method  uses  the  BINTS  calculated  for  the  taped  data.  Coincident  pairs  of  r 
and  Z12  did  not  exist  for  the  telemetered  data  of  lESs  B2  and  B5.  Thus,  the  BINTS  in  column  2  for 
lESs  B2  and  B5  are  estimated  using  the  method  described  in  the  text.  lES  B2  was  never  recovered 
so  there  are  no  entries  in  columns  3  and  4  for  it.  Column  4  shows  the  BINTs  calculated  for  the  data 
put  to  tape.  The  Bints  labeled  with  a  t  indicate  that  their  calculation  by  the  “direct  method”  used 
inferred  rs  .  See  the  text  for  explanation. 


lES 

BINT  1 

Direct 

(m) 

BINT  2 
Comparison 
(m) 

BINT  Tape 
(m) 

Drift  Rate 
“Best  Fit” 
ds/day 

B2 

3367^ 

N/A 

N/A 

+2.16 

B3 

3162 

3229 

7533 

+2.11 

B4 

3457 

3469 

1794 

-1.30 

B5 

3466^ 

3472 

1808 

-2.19 

5.2.3  Despike 


The  despike  program  removes  “bad”  data  values  if  they  fail  either  of  two  criteria.  In  the  first 
test,  data  values  which  fall  outside  some  specified  upper  and  lower  limits  are  “windowed  out”  as 
“bad”  values.  The  range  of  the  “good  data”  window  differed  slightly  for  each  lES  dependent  upon 
the  average  location  of  the  lES  relative  to  the  Gulf  Stream  North  Wall.  For  the  four  lESs,  the 
windowing  range  averaged  1145  ds  (corresponds  to  444  m  main  thermocline  depth).  In  the  second 
test,  data  which  differ  from  a  local  mean  (10  day  running  average)  by  more  than  a  specified  amount 
is  labeled  as  a  “bad”  value.  The  maximum  acceptable  change  of  r  with  respect  to  the  local  mean 
was  set  at  250  ds  which  corresponds  to  a  change  in  Z12  of  97  m. 
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5.2.4  Seacor  (Seasonal  Thermocline  Correction) 

The  program  SEACOR  corrects  the  rs  for  variations  due  to  changes  in  sound  speed  in  the  top  200 
m  of  the  water  column  as  the  seasonal  thermocline  warms  and  cools.  The  seasonal  correction  is  fully 
described  in  Fields  et.  al.,  (1991).  Briefly,  however,  the  rs  are  corrected  using  a  long-term  mean 
curve  of  travel  time  versus  season  calculated  from  over  12,000  XBT  casts  (obtained  from  the  NODC 
archives)  in  the  geographic  region  of  the  lES  deployments.  The  correction  is  important  since  the 
seasonal  variation  of  r  in  the  SYNOP  array  region  was  determined  to  be  roughly  1-1.8  ms  peak  to 
peak  which  corresponds  to  an  apparent  20-36  m  peak  to  peak  variation  in  Z\2  . 

5.2.5  Zi2 

The  Zi2  program  convprts  the  DD.A.CT  directly  into  Z12  depths  using  the  individual  Bints  for  each 
lES  in  Equation  9.  Prior  to  recovery,  the  Bints  in  column  1  of  Table  4  were  used  in  the  program. 
After  recovery  (and  for  all  subsequent  results  shown  in  this  paper),  the  more  accurate  Bints  in 
column  2  for  lESs  B3-B.j  were  used. 

5.2.6  Fill-in 

Gaps  in  the  telemetry  data  records  arise  from  2  sources,  downtime  of  the  Bermuda  receiving  station 
and  false/missed  rs  .  Two  different  methods  were  tested  to  fill  these  gaps  in  order  to  obtain  complete 
continuous  time-series  records.  The  first  method  was  simply  to  linearly  interpolate  between  the  good 
data.  The  results  were  not  very  satisfactory  because  of  the  relatively  large  data  gaps  (—  15  days) 
during  certain  periods.  When  the  data  gaps  were  large,  the  linear  interpolation  scheme  resulted  in 
false  structure  to  the  Z12  record.  The  second  method  proved  a  bit  more  effective  in  filling  the  data 
gaps.  That  method  used  the  spectral  information  present  in  the  data  by  least-squares  fitting  the 
data  to  a  trigonometric  series.  For  example,  let  n  be  the  actual  number  of  data  points  and  .V  be 
the  number  of  data  points  if  there  were  no  missing  data.  The  number  of  missing  data  points  is  then 
m  =  N  —  n.  The  incremental  frequency  is  taken  to  be  (iVAt)~*,  where  At  is  the  sampling  time. 
Hereafter,  At  will  be  taken  to  be  equal  to  1  day.  The  maximum  frequency  is  N(2NAt)~^  =  0  5. 
Letting  h{t)  represent  the  demeaned  data  and  h(t)  =  cos(2x/,t)  fl,  sin(2T/,t)},  th 

method  minimizes  S  =  respect  to  the  Aj  and  Bi.  Since  there  are  missing  da- 

points,  this  system  is  underdetermined  and  solved  by  singular  value  decomposition  (SVD)  and  only 
(Ai,  Bi)  values  are  non-zero  for  non-singular  eigenvalues  of  the  least-squares  fit.  The  fitted  function 
is  then  “back-transformed”  to  yield  the  time  series.  When  large  data  gaps  are  present,  the  method 
fills  them  with  the  (local)  mean  value  of  the  record.  In  practice,  the  despiked  Z12  record  for  each 
lES  was  separated  into  four  approximately  equal  records.  Each  segment  was  run  through  the  fill-in 
program  and  then  the  segments  were  recombined.  Gaps  within  each  segment  were  filled  with  the 


i 
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mean  Zi]  for  that  segment.  The  filled  records  generated  by  the  second  method  are  shown  in  the 
fifth  panel  of  Figures  10 —  13.  The  plateaus  result  from  the  data  gaps. 


5.2.7  Low-pass  Filter 

The  filled  Z12  records  for  the  telemetered  data  are  low-pass  filtered,  with  a  cutoff  period  of  four 
days  to  remove  high  frequency  fluctuations.  The  filter  used  is  a  2nd-order  Butterworth  filter  applied 
forward  and  backwards  to  avoid  phase  shifts.  Note  that  the  taped  data  is  low-passed  filtered  wuh 
a  cutoff  period  of  2  days 


6  Results 

6.1  Data  Return 

On  average,  the  receiving  station  on  Bermuda  was  down  25%  of  the  time.  For  a  given  lES.  the 
downtime  is  calculated  as  any  day  during  which  no  telemetered  data  was  recovered  due  to  problems 
at  the  receiving  station.  This  facility,  however,  was  originally  intended  to  serve  as  a  monitoring 
station  for  the  SOFAR  sound  sources  and  the  downtime  did  not  adversely  effect  that  task  While 
our  data  return  is  much  lower  than  anticipated  because  of  the  downtime,  the  receiving  station  was 
adequate  to  test  the  telemetry  system.  Figure  15  indicates  the  data  return  for  each  lES.  This  data 
includes  that  which  was  recovered  from  the  hard  copies  printed  by  the  ticker-tape  printer. 
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B4;  74% 

B3:  72% 

B2:  78% 

j  1 1  ri'1  TTi  I  [  1 1  I  M  I  1  1  I  1 1 ' :  1 1 1  I  iTj  1 1 — rTTm-p-i — 1 1  ,  1 1  i  j  1  1  1  1 1  1 1  1 1  1  1  m  ■  :  :  m  :  ; 

200.  300  400,  500 

Yearday  ( 1 989) 


iiiiaiaMtit  '  ii^is  tm  iiMia  tiMisiMi*  ■!«  ittfsi 


Figure  15:  The  days  in  which  data  was  received  and  recorded  at  Bermuda  and  received  at  I'RI 


6.2  Dlita  Acceptance 


The  data  acceptance  ia  defined  aa  the  amount  of  data  which  paaaea  the  deapiking  testa  relative  to  the 
received  data.  The  average  data  acceptance  rate  was  79%.  That  is,  on  average,  79%  of  the  received 
data  passed  our  despiking  criteria  involving  windowing  and  maximum  rate  of  change  tests. 

6.3  Delay  Time  Comparison 

The  most  direct  test  of  the  telemetry  system  is  to  compare  the  data  transmitted  by  the  lES  (“data 
sent”)  to  that  received  on  Bermuda  (“data  received”).  The  description  of  how  the  transmitted  data 
is  reconstructed  from  the  taped  data  is  in  Section  5.1.3.  The  variation  in  the  “data  sent"  record 
should  ideally  match  the  variation  in  the  “data  received"  record.  However,  there  will  be  an  offset 
between  the  records  because  the  transit  time  of  the  signal  from  an  lES  to  Bermuda  is  not  precisely 
known.  Another  source  of  offset  is  introduced  because  the  taped  data  truncates  the  5  MSBs  from 
the  original  lES  data  word  while  the  telemetered  d.ita  word  truncates  the  7  MSBs.  For  lES  B3.  the 
offset  can  be  seen  in  Figure  14  to  be  approximately  1 1 150  BDACTs  or  1 1 150  ds.  We  can  remove  this 
offset  by  removing  the  means  from  the  two  records.  Figures  16-  18  show  the  comparisons  between 
the  transmitted  and  received  data  for  sites  B3.  B4.  and  B5.  No  comparison  was  made  for  the  lES 
at  site  B2  since  the  instrument  was  not  recovered.  Both  the  tape  and  the  telemetered  data  shown 
in  these  figures  have  been  despiked. 

Listed  in  Table  5  are  the  root*mean-square  (rms)  differences  between  the  tape  and  telemetered 
records,  listed  both  directly  in  BDACT  and  in  equivalent  meters  of  Zu  depth.  We  used  two  different 
methods  to  compute  these  rms  differences.  In  Method  1,  the  tape  and  telemetered  data  were  first 
independently  despiked  and  the  rms  differences  were  determined  by  subtracting  these  cleaned  up 
records.  This  method  gives  a  true  measure  of  how  well  the  system  works  in  real  time  applications 
because  it  gives  the  rms  difference  between  the  two  data  sets  when  each  is  independently  processed. 
In  Method  2,  we  first  subtracted  the  “raw"  tape  and  telemetry  data  then  despiked  the  difference 
record  before  computing  the  rms.  This  method  gives  values  close  to  Tom  Rossby's  20  ds  scatter  in 
arrival  times  of  the  signals  from  the  moored  SOFAR  float  (see  Section  2.4)  since  the  only  variability 
not  removed  from  the  system  is  the  scatter  due  to  multipath  ambiguities  of  the  transmitted  signal 
between  the  inlet  array  and  Bermuda. 

6.4  Z\2  Comparison 

A  comparison  of  the  Z\%  records  for  all  of  the  telemetry  lESs,  with  the  exception  of  the  one  at 
site  B2,  is  shown  in  Figure  19.  Using  only  days  when  good  data  were  received  at  Bermuda,  the 
rms  difference  between  Z13  as  determined  from  the  tape  data  and  the  telemetry  data  is  29  m  when 
averaged  over  all  of  the  lESs.  Table  6  shows  the  rms  values  for  the  three  individual  lESs,  ranging 
from  23-38  m.  These  values  are  higher  than  those  shown  previously  in  Table  5  where  the  average 
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Figure  17:  Comparison  between  the  time  delay  as  computed  from  the  lES  tape  and  that  measured 
at  Bermuda  for  lES  B4.  Both  records  have  been  independently  despiked  and  demeaned. 
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Figure  18;  Comparison  between  the  time  delay  as  computed  from  the  lES  tape  and  that  measured 
at  Bermuda  for  lES  B5.  Both  records  have  been  independently  despiked  and  demeaned. 


Table  5:  rnu  difierences  of  BDACT  between  the  tape  and  telemetry  data.  The  values  listed  under 
Method  1  are  computed  after  windowing  the  tape  data  and  despiking  the  telemetry  data  separately. 
The  values  listed  under  Method  2  are  computed  after  windowing  out  outliers  in  the  daily  difference 
between  the  t^>e  and  telemetry  data  delay  times. 


lES 

Method  1 

rms 

(ds)  (m) 

Method  2 

rms 

(ds)  (m) 

B3 

40 

16 

r22~ 

9 

B4 

35 

14 

22 

9 

B5 

73 

28 

30 

12 

rms  difference  was  19  m.  These  differences  can  be  attributed  to  the  different  type  of  processing 
performed  on  the  various  data  sets. 

In  the  previous  section,  we  performed  a  “consistency  check"  by  comparing  the  transmitted  data 
with  the  received  data  (Comparison  1 )  We  produced  the  time  series  records  of  transmitted  data  from 
the  taped  data  using  the  same  algorithm  that  the  lES  uses  to  calculate  one  r  for  each  day  and  encode 
it  for  transmission.  Ideally,  the  transmitted  data  and  received  data  should  agree  perfectly.  However, 
noise  is  introduced  when  data  are  transmitted  over  long  distances.  Thus,  the  rms  differences  reported 
in  Table  5  represent  how  well  the  telemetry  system  itself  worked. 

In  this  section,  the  comparison  is  made  between  the  tape  and  telemetered  data  which  have 
been  processed  into  Z13  independently  of  one  another  (Comparison  2).  The  tape  travel  times  were 
processed  using  the  steps  described  in  Section  5.1,  whereas  the  telemetry  data  were  processed  as 
described  in  Section  5.2.  Although  the  processing  steps  perform  essentially  the  same  steps  on  the 
telemetered  data,  they  are  quite  different  in  their  implementation.  For  example,  the  half-hourly  tape 
records  are  detided,  low-pass  filtered,  and  subsampled  to  produce  Zu  time  series  with  a  6-hr  sample 
rate.  Subsequently,  the  tape  Z13  data  are  averaged  to  produce  daily  values.  On  the  other  hand,  the 
data  telemetered  to  Bermuda  is  a  daily  average  and  only  requires  low-pass  filtering.  (Recall  that  the 
taped  and  telemetered  data  are  low-passed  filtered  with  cutoff  periods  of  2-d  and  4-d,  respectively). 
The  ordering  of  the  averaging  and  the  filtering  steps  results  in  the  observed  differences  (Figure  19) 
between  the  tape  and  telemetered  Z13  records. 

It  is  not  unreasonable  to  expect  different  processing  methods  to  yield  different  results.  So  the 
rms  values  in  Table  6  can  be  expected  to  be  higher  than  those  in  Table  5,  where  the  processing  steps 
are  identical.  The  importance  of  the  values  in  Table  6  and  the  records  in  Figure  19  is  that  they 
describe  how  accurately  the  telemetry  system  works  for  monitoring  the  thermocline  depth.  Given 
the  wide  ranges  in  thermocline  depth  at  each  of  the  sites.  Figure  19  shows  that  the  telemetry  system 
works  very  well. 


Depth  of  The  12“C  Isotherm  (m) 


Table  6:  rma  differences  of  Zu  between  the  tape  and  telemetry  data.  Column  one  shows  the  rms 
differences  and  column  two  shows  the  equivalent  differences  in  ds  or  BDACT. 


lES 

Difference  (rms) 

(m)  (ds) 

B3 

38 

98 

B4 

23 

59 

B5 

26 

67 

6.5  Gulf  Stream  Position 


Next  we  used  the  tape  data  and  telemetered  data  to  obtain  the  Gulf  Stream  North  VV’all  (GSNW) 
position  by  2  different  techniques. 

For  the  tape  data,  we  performed  our  standard  method  by  objectively  mapping  Zij  from  the  lESs 
in  the  SYNOP  Inlet  Array  The  GSNW  position  along  a  line  through  the  “B  sites”  was  calculated 
from  these  maps  as  the  position  where  Zn  was  350  m  (note  that  our  standard  definition  of  the 
GSNW  is  where  Zn  is  400  m,  but  Z12  at  350  m  suffices  for  a  comparison  of  the  two  data  sets). 

For  the  telemetered  data  a  more  simplistic  approach  for  determining  the  GSNW  position  was 
used.  Since  the  working  telemetry  lESs  were  moored  along  a  line  approximately  perpendicular  to 
the  Gulf  Stream  path,  a  cross  stream  profile  was  used  to  determine  the  GSNW  position  using  the 
technique  of  Watts  and  Johns  (1982). 

A  comparison  between  the  Gulf  Stream  North  Wall  (GSNW)  position  calculated  from  a  common 
geographic  position  for  both  the  tape  and  telemetered  data  is  shown  in  Figure  20.  For  the  telemetered 
data,  on  any  given  day,  the  GSNW  position  is  determined  from  some  combination  of  actual  and  filled 
Z12  values.  The  open  triangles  in  Figure  20  show  the  GSNW  position  on  days  when  there  was  at 
least  one  lES  reporting  “good  data”  (in  this  case  the  GSNW  position  is  calculated  from  1  actually 
measured  Zu  and  3  filled  Zi2  va;values).  Conversely,  the  solid  triangles  denote  days  when  the  GSNW 
position  was  calculated  from  4  filled  Z12  values.  The  low  frequency  signal  shows  good  agreement 
except  during  the  periods  when  the  receiving  station  was  down  while  the  Gulf  Stream  moved  rapidly 
onshore  (—  yeardays  255-265  and  267-275).  The  rms  difference  is  less  than  5  km  between  the  tape 
and  telemetered  records  for  the  Gulf  Stream  position  for  days  in  which  good  data  was  received 
from  at  least  one  lES.  The  5  km  rms  is  better  than  the  ~10  km  GSNW  resolution  from  SST  maps 
generated  from  AVHRR  data.  Thus,  the  telemetry  system  is  an  excellent  method  of  obtaining  real¬ 
time  GSNW  position  through  a  limited  domain  mapped  with  lESs  measuring  main  thermocline 


6.6  Telemetry  System  Problems 

This  WM  the  first  deployment  of  a  telemetry  system  for  real-time  aquisition  of  lES  data  and  the 
system  concept  was  proven  to  work.  However,  as  is  frequently  the  case  when  things  are  done  for  the 
first  time,  a  few  problems  occurred  for  which  we  have  some  proposed  solutions. 

As  detailed  earlier  in  Section  5.2.2,  a  very  significant  problem  with  the  telemetry  system  occured 
with  the  frequency  settings  of  the  lES  internal  clocks.  Our  frequency  counter,  used  to  measure  the 
clock  frequencies,  was  not  calibrated  correctly  by  the  factory  so  our  frequency  measurements  were 
in  error.  As  a  result,  the  telemetry  records  had  time  drifts  which  were  translated  into  Zu  drifts  by 
the  data  encoding  scheme.  The  data  could  not  be  satisfactorily  de-drifted  until  the  lES  data  tapes 
were  recovered.  This  drift  problem  could  have  been  avoided  by  including  a  periodic  zero-delayed 
signal  from  each  lES.  This  signal  could  be  used  to  keep  track  of  any  clock  drift  since  such  a  drift 
would  manifest  itself  as  a  departure  of  the  arrival  time  of  the  signal  from  the  expected  time. 

Gaps  in  the  data  were  the  result  of  both  received  data  which  were  “bad"  (rejected  by  the  tesi.s 
for  bad  data  mentioned  in  section  .5.2.3)  and  receiving  station  downtime.  From  the  top  panels  of 
Figures  10-13  it  is  seen  that  the  “bad”  data  were  concentrated  during  a  70  day  time  span  centered 
about  yearday  360  of  1989.  That  is,  the  noisiest  data  were  during  the  months  of  December  through 
February,  the  time  period  of  severest  weather  over  the  Sargasso  Sea  and  Bermuda.  At  present,  the 
only  remedy  for  this  problem  is  to  interpolate  between  good  values  to  fill  in  the  gaps  caused  by 
the  rejection  of  these  points.  The  gaps  resulting  from  the  receiving  station  downtime  are  something 
that  could  have  been  almost  entirely  eliminated.  Until  an  uninterruptable  power  supply  (UPS)  was 
added  to  the  system,  data  were  lost  whenever  the  power  went  down  (the  power  supply  on  Bermuda 
is  not  very  stable).  Data  were  also  lost  when  problems  arose  with  the  modem  used  in  the  system. 
As  pointed  out  in  Section  6.1,  the  downtime  of  the  receiving  station  was  about  25  %  so  there  is 
significant  room  for  improvement. 
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Figure  20:  Gulf  Stream  North  Wall  displacement  calculated  from  Objectively  Mapped  fields  of  Zi2 
in  the  SYNOP  inlet  array  along  the  “B  line”  of  lESs  and  from  the  telemetered  data.  The  calculations 
are  referenced  to  a  common  geographic  position. 
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7  Summary  and  Conclusions 

Despite  the  few  avoidable  problems  encountered  in  the  system,  it  has  been  shown  that  it  is  possible 
to  obtain  high  quality  lES  data  using  an  acoustic  telemetry  system.  The  average  rms  difference  of 
19  m  between  the  data  telemetered  from  the  lESs  and  the  data  received  at  Bermuda  is  nearly  the 
same  as  the  average  20  m  uncertainty  in  BINT  calculations  for  lESs  in  our  study  area.  Fields  et 
al.,  (1991).  This  system  was  used  for  real  time  tracking  of  the  Gulf  Stream  with  accuracy  of  o  km 
The  telemetry  system  is  not  limited  to  lESs,  but  could  include  similarly  processed  information  from 
current  meters,  pressure  gauges,  ambient  noise,  etc.  Such  a  system  could  be  used  in  other  regions 
where  lESs  are  used  (such  as  the  Kuroshio)  where  there  is  a  well  defined  sound  channel 
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